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Abstract 

We study parameter estimation with post-Newtonian (PN) gravitational waveforms for the quasi- 
circular, adiabatic inspiral of spinning binary compact objects. In particular, the performance 
of amplitude-corrected waveforms is compared with that of the more commonly used restricted 
waveforms, in Advanced LIGO and EGO. With restricted waveforms, the properties of the source 
can only be extracted from the phasing. In the case of amplitude-corrected waveforms, the spectrum 
encodes a wealth of additional information, which leads to dramatic improvements in parameter 
estimation. At distances of ~ 100 Mpc, the full PN waveforms allow for high-accuracy parameter 
extraction for total mass up to several hundred solar masses, while with the restricted ones the 
errors are steep functions of mass, and accurate parameter estimation is only possible for relatively 
light stellar mass binaries. At the low-mass end, the inclusion of amplitude corrections reduces the 
error on the time of coalescence by an order of magnitude in Advanced LIGO and a factor of 5 
in EGO compared to the restricted waveforms; at higher masses these differences are much larger. 
The individual component masses, which are very poorly determined with restricted waveforms, 
become measurable with high accuracy if amplitude-corrected waveforms are used, with errors as 
low as a few percent in Advanced LIGO and a few tenths of a percent in EGO. The usual spin-orbit 
parameter jS is also poorly determined with restricted waveforms (except for low-mass systems in 
EGO), but the full waveforms give errors that are small compared to the largest possible value 
consistent with the Kerr bound. This suggests a way of finding out if one or both of the component 
objects violate this bound. On the other hand, we find that the spin-spin parameter a remains 
poorly determined even when the full waveform is used. Generally, all errors have but a weak 
dependence on the magnitudes and orientations of the spins. We also briefly discuss the effect of 
amplitude corrections on parameter estimation in Initial LIGO. 
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I. INTRODUCTION 



Inspiraling compact binary objects (neutron stars and/or black holes) are among the most 
promising sources for existing and planned interferometric gravitational wave observatories 
[l|]. Efforts are in progress to detect their gravitational wave signals in LIGO, VIRGO, 
GEO600 and TAMA [2| . When the gravitational wave signal enters the bandwidth of these 
detectors, typically the orbits of the bodies will have approximately circularized, and the 
subsequent inspiral is described well by an adiabatic approximation in which the radii of the 
orbits are assumed not to change significantly over a single period [3|. In this regime the 
frequency and amplitude of the signal increase steadily in a "chirp" until a last stable orbit 
(LSO) is reached. The amplitude and phasing of signals arising from quasi- circular adiabatic 
inspiral have been computed to great accuracy using post-Newtonian (PN) methods (see 
for a review), which lead to expressions that take the form of series expansions in the orbital 
velocity v. 

Apart from their obvious astrophysical importance, binary inspirals are potential "lab- 



oratories" for testing general relativity This is because at least in the adiabatic 

phase they are relatively "clean" and reasonably well-understood systems. However, it is 
important to find out how good a testing ground they offer. Precisely what information do 
the waveforms carry? Which parameters associated with the binary can be determined from 
them, and to what accuracy can this be achieved? In other words, to what extent will it be 
possible to reconstruct binary inspiral events from the gravitational radiation they emit? 

The most accurate waveforms currently available for the quasi-circular, adiabatic inspiral 
of compact objects are of order 0{v^) in amplitude [sl and C(w^) in phase [9|, at least in the 
case of non-spinning binaries; in the usual notation this corresponds to 2.5PN and 3.5PN 
orders, respectively. An efficient way of searching for inspiral signals in data is matched 
filtering [10|, which involves a bank of templates. Post-Newtonian waveforms describing 
adiabatic inspiral are linear combinations of harmonics in the orbital phase. In most of 
the literature, simulated signals as well as detection templates are taken to be restricted 
post-Newtonian waveforms , which only have the dominant harmonic at twice the orbital 
phase, with PN corrections to the phasing but no corrections to the amplitude. The effects on 
signal-to-noise ratio (SNR) and redshift reach as a result of including amplitude corrections 
in templates and simulating signals have been investigated in 1^, 1^ in the context of both 



initial and advanced detectors, and it was found that the consequences are considerable. In 
initial detectors the use of restricted PN waveforms can lead to a significant overestimation 
of SNR. More interestingly, in advanced detectors the higher harmonics can give rise to 
considerable SNRs even when the dominant harmonic (and hence the restricted waveform) 
does not even enter the detector's bandwidth. As a result, the mass range of such detectors 
increases dramatically with the inclusion of amplitude corrections in template waveforms. 
Also in the case of masses for which the dominant harmonic reaches the bandwidth in only 
a small frequency interval, taking into account all known harmonics greatly increases the 
SNR. Thus, for high stellar mass to intermediate mass binaries (100 Mq ^ M < 1000 Mq), 
detection rates increase dramatically when amplitude corrections are taken into account. 
This also opens up the possibility of high-quality parameter estimation up to much higher 
masses than would be thought possible on the basis of the restricted PN approximation. 

More generally, to properly address the questions raised above, we will see that it is 
crucial to use the amplitude-corrected waveforms to highest available post-Newtonian order. 
With restricted waveforms /iq, almost any information one can extract about the source has 
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to come from the phasing, i.e. arg(/io(/)), with Hq the Fourier transform of the time domain 
waveform. In the case of amphtude-corrected waveforms h, on top of the information carried 
by arg(/i(/)) one can exploit the intricate structure of the spectrum, i.e. the power per 
frequency bin as a function of frequency. As we shall see, the additional information carried 
by the spectrum of the full waveform has dramatic consequences for parameter estimation, 
whence the title of this paper. (The study of black hole ringdown in terms of quasi-normal 



modes is often also referred to as "spectroscopy" [ij]; in the case of inspiral the spectrum 
is of course continuous.) 

Most of the research to date on parameter estimation with inclusion of additional har- 
monics has been of a preliminary nature. For Initial LIGO there is the work of Sintes and 
Vecchio [3], who used waveforms at 0.5PN order in amplitude and 2PN order in phase. 
Very recently. Rover et al. [3] considered a network of initial detectors (the Initial LIGO 
detectors and VIRGO), taking the sources to be binary neutron star inspirals and using 
waveforms at 2PN in amplitude and 2.5PN in phase. For the proposed space-based detector 
LISA there have been studies by Sintes and Vecchio [13] (again at 0.5PN in amplitude and 
2PN in phase) and by Hellings and Moore [3, [l^ (2PN in both amplitude and phase). 
None of these works took the effects of spin into account. In this paper we focus on ad- 
vanced ground-based interferometric gravitational wave detectors. We will mostly consider 
parameter estimation with Advanced LIGO as well as a European third-generation detector 
that is currently under consideration [13]; the latter will tentatively be called the Euro- 
pean Gravitational- wave Observatory (EGO). The waveforms used will be of 2.5PN order in 
amplitude and phase, with inclusion of spin-related contributions up to 2PN in the phase. 

Up to 2.5PN order in amplitude, the waveforms in the two polarizations take the general 
form [sl 

h^,. = {<i + xV^</f) + xH^^l + x'/'H^^;^^ + x'Hf^, + } (1.1) 

where r is the distance to the binary, M its total mass, and r] the ratio of reduced mass 
to total mass. We have set G = c = 1, as we will do throughout this paper. The post- 
Newtonian expansion parameter is defined as x = (27rMF)^/^, with F(t) the instantaneous 

orbital frequency. The coefficients H^^^\ p = 0, . . . , 5, are linear combinations of various 
harmonics with prefactors that depend on the inclination angle t of the angular momentum 
of the binary with respect to the line of sight as well as on rj] their explicit expressions can 
be found in [81]. The measured signal also depends on the polarization angle and the position 
in the sky through the detector's beam pattern functions F^^x- 

h{t) = F^h+{t) + Fxhx{t). (1.2) 

Note that for ground-based detectors, which are the ones we will be concerned with, it is 
reasonable to approximate -F+,x as being constant over the duration of the signal. They 
depend on angles {6,(f),ip), where (6*, 0) determine sky position while is the polarization 
angle. The signal (11.21) is a linear combination of harmonics of the orbital phase \l/(t) with 

offsets ip(k,m/2y- 

N 5 

h{t) = X] Ak,m/2){t) C0s(/C^'(t) + ^P(k,m/2)), (1-3) 

fc=l m=0 

where k runs over harmonics while m/2 is PN order in amplitude. The coefficients A(^k,m./2) 
are functions of (r, M, r], 9, 0, ip, l) multiplied by x*^'""'"^)/^, with x depending on time through 
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F{t). The number of harmonics depends on the PN order in amphtude; at 2.5PN one 
has N = 7. 

Spin-related effects will be taken into account in the following way. In the phasing, spin- 
orbit effects are included at 1.5PN and 2PN orders, and spin-spin effects at 2PN order, as 
in j2l|. (The phasing formula with the inclusion of spin-orbit interactions up to 2.5PN has 
recently been computed [13], but the research we report on here was already well underway 
when these results appeared.) We will use the stationary phase approximation to the Fourier 



transform [23|], in which the coefficients of the harmonics are proportional to the inverse 
square root of the time derivative of the associated instantaneous frequency. We let spin- 
related effects enter the amplitudes through these "frequency sweeps" , at the same respective 
orders. Since we want to view the full waveform as a model for an actual signal as well as 
a search template, for consistency we take both amplitude and phase to be of 2.5PN order. 
This is a choice we make, motivated by conceptual considerations; there are no technical 
problems in taking the phasing up to 3.5PN. 

Both the expressions for the amplitudes in (11.31) and the expression for the phase can 
be truncated to a lower post-Newtonian order. Such truncated waveforms will be denoted 
as [p, q)PN, where p refers to the PN order of the amplitude and q that of the phase. 
Thus, the restricted PN waveform will be referred to as (0, 2.5)PN and the full waveform as 
(2.5,2.5)PN. 



For the Initial LIGO detector, Sintes and Vecchio [15| discussed the improvement in error 



estimation gained by using amplitude-corrected waveforms by means of the covariance matrix 



formalism [2J]; in the present paper the same techniques will be employed. Concretely, 
the following questions will be addressed: (i) Using the best available waveforms, which 
parameters can be determined? (ii) How do the estimation accuracies improve as the PN 
order of the amplitude is increased? and (iii) What is the effect of the magnitudes and 
orientations of the components' spins?^ 

This paper is structured as follows. In section [III we first briefly review the covariance ma- 
trix formalism for parameter estimation. The Fourier transform of the amplitude-corrected 
waveforms is computed in the stationary phase approximation and its spectrum is discussed. 
We then apply the covariance matrix formalism to (p, 2.5)PN waveforms with p > 0.5. In 
section IIIII we compare the performance of restricted and amplitude-corrected waveforms 
in EGO and Advanced LIGO for a variety of systems. We also study the improvement in 
error estimation by using waveforms of order [p, 2.5)PN as p varies from zero to 2.5, and 
the dependence of errors on the spins of the component objects. Finally, we present results 
for some representative systems with Initial LIGO. Section IIVI provides a summary and 
conclusions. 



^ Throughout this paper we will include a spin-orbit and a spin-spin parameter as variables in the error 
analysis so that even when spins are set to zero, the corresponding errors are never neglected. 
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II. AMPLITUDE-CORRECTED POST-NEWTONIAN WAVEFORMS AND PA- 
RAMETER ESTIMATION 



A. Overview of the covariance matrix formalism 



We begin by recalling some basic facts of the covariance matrix formalism for parameter 
estimation with matched filtering [1^. To search for signals buried in the noise, the data 
analyst uses families of waveforms h[6"']{f), determined by a finite number of parameters 

that characterize the binary. The measured values of the parameters, O"", depend on 
the realization of the detector noise at the time the signal was present, so that in general 
they do not correspond to the "true" values 6°'. For large SNR, the measurement errors 



A6'" = O"' — 6"' follow a multivariate Gaussian probability distribution [24 



P(Ar) = Po exp(-r6e(^'') A^^A^72) (2.1) 

with Po ^ normalization factor, and Tab{d"') is the Fisher information matrix. The latter 
may be thought of as a metric on the differentiable manifold naturally defined by the family 
of waveforms {/i[6'"]}, and it is given by 

Tab = {ha\h), (2.2) 

where ha = dh/d6°', and ( . | .) is the usual product defined by the one-sided noise power 
spectral density Sh{f)- 

ff-- Re{i*{f)y{f)) 




The integration domain .fmn.r] is determined both by the detector and by the nature 

of the signal. As explained in [13| (and in subsection III Bl below) . it is reasonable to cut off 
the kth. harmonic in the waveform at a frequency kfLso, where /lso may be viewed as the 
orbital frequency at the last stable orbit. Since the amplitude-corrected waveform contains 
7 harmonics, in that case we will effectively have fmax = ^fiso, while for the restricted 
waveform fmax = '^fbso- As to the lower cut-off frequency, power spectral densities Sh{f) 
tend to rise very quickly below a certain frequency fg where they can be considered infinite 
for all practical purposes; accordingly, we take fmin = fs- 

The signal-to- noise ratio is also defined in terms of the product (12. 3p . 

p[h] ^ ^{h\h)- (2.4) 
The inverse of the information matrix is the covariance matrix S'^*, and one has 

^ab ^ ^-p^lyb ^ ^ArA^^) (2.5) 

where (.) denotes the average with respect to the probability distribution P(A6'") (12.11) . 
The root-mean-square error in the measurement of the parameters is then 



a" = V((A^»)2) = (2.6) 

and with slight abuse of notation these will simply be denoted A0° in the rest of the paper. 
Finally, the correlation coefficients are defined by 

c-b ^ (2.7) 



5 



(Note that repeated indices in Eqs. (12. 6p and (12. 7p are not meant to be summed over.) By 
definition one has — 1 < Cat < 1. If the correlation coefficient between two parameters is 
close to 1 (resp. —1), they are strongly correlated (resp. anticorrelated), meaning that one of 
them is redundant. A correlation coefficient being close to zero indicates that the parameters 
are largely unrelated. 

Usually the Fisher information matrix, and hence errors and correlations, can not be 
calculated exactly because the integrals defining Tab can not be performed analytically. This 
is also the case here, and most of the results in this paper will refer to numerical calculations 
with the software package Mathematica. The matrix E = F"^ was obtained using the built- 
in matrix inversion routine. In that regard it is worth mentioning that problems can arise 
when numerically inverting large matrices such as the Fisher matrices encountered here; we 
will come back to this point in subsection IIIDI below. 

The covariance formalism is only valid for sufficiently high SNR; at low SNR it tends 
to underestimate errors, as seen in Monte Carlo treatments involving large numbers of 
simulations with different realizations of the noise [1^. However, at distances of a few 
hundred Mpc the SNRs in Advanced LIGO and EGO will be rather high, as was made 



evident in 1311: see also subsection IIII Al below. 



B. Amplitude-corrected PN waveforms 

The Fisher matrix (12. 2p is defined in terms of (derivatives of) the Fourier transform of 
the waveform h[6"']{t). We will find it convenient to use the well-known stationary phase 
approximation (SPA) [2^ to the Fourier transform in order to compute the scalar products 
involved. The resulting waveforms have already been discussed in detail in l^]; here we 
only give a quick overview. 

As mentioned in the Introduction, the waveforms h[6"-]{t) take the form 



h{t)=J2h^'\t) (2.8) 

k=l 

where the h^''"' (t) involve harmonics of the orbital phase with constant offsets: 

5 

). (2.9) 

m=0 

These harmonics are found from (11.11) and (II. 2p : sky position [6, (p) and the polarization angle 
ip enter the coefficients ^(fc,m/2) and offsets (p(k,m/2) through the beam pattern functions F+,x: 

F+{e,(f),ip) = i(l + cos^(e))cos(20)cos(2V^)-cos(e)sin(20)sin(2V'), 
F^{e,(l),ij) = ^ (1 + cos2(e)) cos(2(/))sin(2V^) + cos(^)sin(20)cos(2V^), (2.10) 

which in the case of ground-based detectors can be considered constant for the duration of 
the observed part of the signal. 

During the inspiral phase one has \d\nA(^k^s-j/dt\ <C 1 and \kd'^'^ /dt'^\ <^ (kd'^/dt)^, in 



which case we can use the stationary phase approximation [23|] to the Fourier transform of 
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fl2.9p . For positive frequencies we have 



E 



m=0 ^(fc,m/2) 



exp 



^ 1 27r/t ( -/ 



q 7/ - <^(fc,m/2) - 7r/4 



exp 



z ( 27r/t, - 7r/4 + k^j [ -f 



(2.11) 



A dot denotes derivation with respect to time and t^. is the coalescence time. 

We recall that the function F{t) is the instantaneous orbital frequency. The function 
t{f) is defined implicitly by F(t(f)) = f. In the expressions (EUD for the h^''^ we take the 
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-r]j {2'KMFf'^ 

(2.12) 

where M. = Myf/^ is the chirp mass. The parameter jS occurring at 1.5PN order encapsulates 
the leading-order spin-orbit effects, while the parameter a at 2PN order also encodes spin- 
spin effects; at 2.5PN order we neglect all spin-related effects. (3 and a are given by [30|]: 



/3 



1 2 

- ^ [113 (mi/M)2 + 75 r/] L ■ x., 



12 

H 

48 



i=l 



-247(xi-X2) + 721(L-xi)(^-X2) 



(2.13) 
(2.14) 



with Xi = Si/mf and Si, i = 1,2, the spins of the binary's components, while L is the unit 
vector in the direction of orbital angular momentum. To a reasonable approximation, (3 and 
a may be considered constant 33|]. The maximum possible value of |(j| for sub-extremal 
black holes is 79/32 ~ 2.5 and the maximum value of \(3\ is 113/12 ~ 9.4. 
To 2.5PN order, the phase il'if) is given by 



(2,15) 



where ijjc is the orbital phase at coalescence. We take the coefficients ijji to be j^, 0, [if 
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The SPA for the {p, 2.5)PN waveform is then 



ISPA 



k=l 



E 



2p 



m=0 ^(fc,m/2) 



g-«¥'{fc,m/2) 



exp 



271 ft, - 7r/4 



(2.17) 

where Np is the number of harmonics, and [ . ]p denotes consistent truncation to pth post- 
Newtonian order (i.e., the "Newtonian" prefactor /"^''^ is taken outside and the remaining 
expression is Taylor-expanded in up to order {2'KMff/^). For more exphcit 

expressions we refer to [l3|. 

In the time domain it is reasonable to cut off the waveform at a time corresponding to 
F{t) = fiso, where /lso is the orbital frequency at "last stable orbit". For simplicity we 
use the expression corresponding to the extreme mass ratio: 



LSO 



63/227rM' 



(2.18) 



where M is the total mass of the system. In the frequency domain this roughly corresponds 
to terminating the harmonics at multiples of Jlso] in practice we multiply the truncated 
kth harmonic by 0{kfLso ~ f) where 6{x) is the usual Heaviside function {6{x) = 1 if x > 
and 9{x) = otherwise) ^ 



Finally, in [28|, |29[ the stationary phase approximation was compared with an alternative 



way of approximating the Fourier transform, the so-called fast Fourier transform (FFT), for 
initial detectors and relatively small masses. In that case the SPA and FFT differ by only a 
few percent, and the situation improves dramatically with increased detector sensitivity at 
low frequencies. We also note that the conditions for the applicability of the SPA as spelled 
out above are satisfied more easily with increasing k. 



C. Spectrum as observed in a detector 

Our goal is to compare the quality of parameter estimation with restricted and amplitude- 
corrected waveforms. As we now explain, the extent to which the full waveform will perform 
better is determined by the information content of its spectrum. Consider the contribution 
to the SNR squared per logarithmic frequency bin, as a function of frequency. The squared 
signal-to-noise ratio of a waveform h m a, detector with noise power spectral density Sh{f) 
is 



mm 



fmii 

We now define the "observed spectrum" as 



Sh{f) 

/IM/)P 
Sh{f) 



d\n{f). (2.19) 



Vif) ^ (2.20) 

The name is apt; the observed spectrum bears a direct relationship to the way signals are seen 
in a detector, which depends both on the waveform and on the sensitivity of the detector. In 
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[29| the quantity V{f) was used (under another name) as a diagnostic to compare different 
(restricted) waveform approximants. 

When using restricted waveforms for parameter estimation, of necessity almost all 
information one can extract about the source has to come from the phasing, arg(/;,o(/)); 
the spectrum does not play much of a role. The situation is very different when amplitude- 
corrected waveforms h are used. In that case we can exploit not only arg(/i(/)); on top of 
that there is the wealth of information contained in |/i(/)p, or more directly in the observed 
spectrum 'P(/). 




200 400 600 800 1000 1200 1400 20 40 60 80 100 120 140 

f (Hz) f (Hz) 



Figure 1: The observed spectrum in Advanced LIGO for a (1,1O)M0 system (left) and a 
(10, 100) Mq system (right). Distance is 100 Mpc and angles were chosen arbitrarily as 6 = (p = ir/d, 
ip = 7r/4, i = tt/3. Spins are zero. 

In Fig.[T]we have plotted V{f) for both restricted and amplitude-corrected waveforms, for 
two different systems in Advanced LIGO. With a (1, 10) Mq system, the restricted waveform 
ho win penetrate a detector's bandwidth up to a relatively high frequency, as 2fiso = 
l/(6^/^7rM) will be large. The higher harmonics that are present in the full waveform 
continue to still higher frequencies, but their contribution beyond 2fLso is suppressed by 
the denominator Sh{f) in f l2.20p : Advanced LIGO has its peak sensitivity at about 215 Hz. 
On the other hand, for a (10, 100) Mq system, ho will deposit power in a much smaller 
frequency interval, with an upper cut-off frequency much below the frequency where the 
detector is most sensitive. Therefore, for the restricted waveform, the quantity V{f) will 
be non-zero in only a small interval in which the denominator Sh{f) in fl2.20p is large. By 
contrast, the higher harmonics in the amplitude-corrected waveform h can penetrate the 
bandwidth up to frequencies where the detector is much more sensitive; note the distinctive 
features at 40 Hz, 60 Hz, 80 Hz, . . .. 

For both systems, V{f) for the full waveform exhibits a great deal of structure due to 
interference between the various harmonics. As a consequence, we may expect the amplitude- 
corrected waveforms to allow for parameter extraction with significantly smaller errors than 
the restricted waveforms, with the difference being much more pronounced for higher masses. 
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D. Parameter estimation with amplitude-corrected waveforms: Overview 



Before going into details, we begin with an overview and discuss parameter estimation 
with the waveforms of the previous subsection in a quahtative way. 

For any family of waveforms of order (p, g)PN with 0.5 < p < oo, the metric Tab has a 
coordinate singularity at 77 = 1/4. This can readily be inferred from the explicit expression 
for (0.5, g)PN waveforms in [l5|: /iin(»?) blows up as 77 1/4 and the same is true of rin(r,)fe = 
(/iin(r,)|/ife). The problem does not occur if the usual coordinate ln(?7) is replaced by 6, where 



6 = 



nil — ^2 



mi + m2 



(2.21) 



(This new parameter was already introduced in [18|, but there the motivation was left 
unclear.) For this reason, when considering amplitude-corrected waveforms we will use the 
set of parameters 

(r) = (ln(l/r),ln(^),(5,t„^e,cos(O,cos(^),0,^,/3,cT). (2.22) 

Generally, not all parameters affecting the waveform can be extracted a posteriori, and 
this will also be the case here. As the PN order in amplitude is increased, one has the 
following: 

• The restricted [p = 0) waveform, whose SPA takes the form 

Ar'/' exp [z (27r/t, - n/4 + 2^/;(//2))] , (2.23) 
only allows for the determination of [1^, 3ll 



{ln{A),ln{M),6,t,,^,,/3,a), (2.24) 
where the amplitude A depends on r, A4, 6, and the four angular parameters. 

Going to 0.5PN in amplitude, two additional harmonics appear and it becomes possi- 
ble in principle to discriminate between three variables appearing in the amplitudes: 
The distance r, the inclination angle l, and e.g. the quotient F^/F^ of the beam 



pattern functions [15[. Below we will instead include (ln(l/r), cos(t), cos(6')), as these 



parameters have a direct physical meaning. 

At IPN in amplitude there are four harmonics in all so that one might expect either 
or ip to be an additional measurable parameter. However, whichever one chooses, the 
errors on (cos(^), 0, ijj) tend to be extremely large and any of these parameters will be 
unmeasurable in practice. The same goes for the distance r. This remains true up to 
2.5PN order in amplitude. 



Ideally one would like to work with the full 11 x 11 Fisher matrix based on all of the 
variables in fl2.22|] . However, inversion of such a large matrix tends to fail when using the 
built-in inversion routine of Mathematica. The reason is that the matrix will often be "ill- 
conditioned" ; this problem was also encountered in |5|, and a good discussion can be found 
in Appendix B of that paper. We were able to avoid this problem by instead considering a 
Fisher matrix that takes into account a smaller number of variables but leads to accurate 
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estimates for the errors on the parameters one would be the most interested in. As discussed 
above, all variables in (12.221) except for r and (cos(6'), 0, tp) can be determined at the highest 
post-Newtonian order considered in this paper. As we shall see, the latter three are strongly 
(anti-)correlated among each other and have only tiny correlations with all other parameters. 
On the basis of this one might conclude that it is safe to disregard them entirely and build 
an 8 X 8 Fisher matrix associated with the eight remaining parameters, 

(ln(l/r), ln(A<), 6, t„ ^j,, cos(0, /?, a). (2.25) 

By comparing results for a representative sample of points in parameter space where inver- 
sion of both the 8x8 and 11x11 Fisher matrices did succeed, we were able to ascertain 
that it is indeed safe to use the smaller matrix, except for the following important caveat. 
Use of the 8x8 Fisher matrix leads to only a few percent difference in the errors on 

(ln(Al),5,te,cos(0,At^), (2.26) 

but it leads to a deceptively small error on distance (Ar/r ~ 0.1 for a large range of masses 
in Advanced LIGO and about ten times smaller in EGO). The reason is that excluding 
{cos{9), (f),ip) from the analysis is equivalent to assuming that they are known, fixed pa- 
rameters. Since is an overall prefactor in the waveform and because A4 is already 
well-determined from the phase alone (26l . 3l| . the error on distance then also comes out to 



be small. In reality, with inclusion in the analysis of sky position and polarization angle, 
the relative error on distance tends to be far larger than 100%. This is as one would expect 
for short-lived sources in a single interferometer j32|. Unless stated otherwise we will work 
with the 8x8 Fisher matrix corresponding to the eight variables (12.251) . which leads to a 
robust estimation of the parameters (12.261) . and the values for the errors on ln(l/r) and ijjc 
will be disregarded. For most choices of parameter values, this smaller Fisher matrix is free 
of the abovementioned problem with matrix inversion. 

For restricted waveforms, the 6x6 Fisher matrix based on the last six parameters^ in 
(I2.24P becomes ill-conditioned when mi/m2 — 1. The reason is that, unlike with the full 
waveforms, with restricted waveforms the error A(5 grows very quickly as mxjm^ ^ 1. To 
understand how this comes about, consider a coordinate system in which ln(?7) is used in 



place of as in |26|. One has A?]/// = \d\a{r]) j db\ , and \d\\\{r]) j db\ = 25/(1 — 5^) ^ 
as mi/m2 — >■ 1, while Arj/r] remains finite; hence AS — >■ oo. The result is that in our 
framework, we will not be able to compare parameter estimation with restricted versus 
amplitude-corrected waveforms when mi/m2 — 1- Indeed, in that case the Fisher matrix for 
restricted waveforms is ill-conditioned if S is used as a coordinate, while the Fisher matrix 
for the full waveforms becomes singular at mi = m2 if ln(?7) is used in place of 6. 

With amplitude-corrected waveforms, the mass-related parameters Ai and 6 and also 
the coalescence time tc can be measured with great accuracy for a wide range of systems. 
One tends to have Acos(t) ~ 0.2 with Advanced LIGO and ten times smaller with EGO. 
However, in the absence of information on the polarization angle the information carried by 
the inclination angle is not very useful, and we will not put much emphasis on it. 

Finally, a comment on the spin- related parameters /? and a. In previous works [sH, 13 
one would multiply the error distribution (12.11) by Gaussians in /5 and a with appropriate 



^ The Fisher matrix based on all of the parameters (|2.24p will be block-diagonal, and the block associated 
with the variable In(^) can be ignored 
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widths to reflect the fact that the absolute values of these parameters have very specific 
upper bounds if the components of the binary are sub-extremal Kerr black holes. In the 



context of restricted PN waveforms it was shown in [31[ that doing so has the effect of 
decreasing error estimates for {\n{Ai),ln{ri),tc,'ipc)- However, apart from the fact that the 
method is somewhat inelegant, it also inculcates our biases as to what the endpoint of 
gravitational collapse should be. Gravitational wave astronomy will provide us with an 
opportunity to directly search for objects whose spins exceed the Kerr bound, as we will 
discuss below. Accordingly we will not put any a priori constraints on the values of (3 and a. 
The parameter f] is well determined in both Advanced LIGO and EGO; a is unmeasurable 
in Advanced LIGO and it is rather poorly determined also in EGO. 

In summary, even with (2.5, 2.5)PN waveforms the number of parameters that can be 
determined with reasonable accuracy is the same as for the restricted ones. With restricted 
waveforms the well-determined variables are (In(^), In(A^), S, tc, (3, ex); with the (2.5, 2.5)PN 
waveforms they are (\n{Ai), 6, cos(t), /3, a). 



III. RECONSTRUCTION OF BINARY INSPIRAL EVENTS USING ADVANCED 
DETECTORS 

We now study the quality of parameter estimation to be expected with Advanced LIGO 
and EGO, comparing the performance of restricted and amplitude-corrected waveforms for 
two sets of binary systems. In one set, the mass ratio qm = mi/m2 is held fixed while the 
total mass M is varied; in the other, M is kept the same but qm is varied. The performance 
of restricted PN waveforms will be compared with the more general, (2.5, 2.5)PN ones. We 
will mostly be interested in asymmetric systems where one component is much heavier than 
the other since it is for such systems that we expect the higher-order amplitude corrections 
to make the biggest difference. When varying M for fixed q^,, for the most part we restrict 
attention to systems with total mass up to ~ 250 Mq in Advanced LIGO and ~ 500 Mq in 
EGO. (Moderately heavy and moderately asymmetric systems with masses in the order of 
hundreds of solar masses are likely to be formed in galactic nuclei in the process that leads 
to the formation of supermassive black holes, or in globular clusters {13].) We will also look 
at the dependence of errors on the PN order in amplitude. Initially the spins of the binary's 
components will be set to zero so that f3 = a = 0; near the end of this section we investigate 
how errors depend on the values of the spins. 

The set of parameters (12.221) can naturally be divided into two subsets: 

(i) The "extrinsic" parameters (tc, "^c, ln(l/r), cos(6'), 0, cos(i), ■?/') which describe 
the time at which coalescence occurred, the phase at coalescence, the distance, 
the sky position, and the orientation of the inspiraling binary with respect to 
the detector; 

(ii) The "intrinsic" parameters {ln[Ai),6, /3,a) which encode the masses and 
spins of the binary's components. 

We discuss in turn the behavior of signal-to-noise ratios, errors on extrinsic parameters, 
and errors on intrinsic parameters. 
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Figure 2: Plots of p[fi\ and p[/io] in Advanced LIGO and EGO as functions of total mass with 
fixed qm = mi/m2 = 0.1 (left), and as functions of Qm with fixed total mass M = 100 Mq 
(right). Distance is fixed at 100 Mpc and spins are set to zero. Angles were chosen arbitrarily as 

e = (/) = 7r/6, ^ = 7r/4, i = vr/S. 



A. Signal-to-noise ratios 



The issue of signal-to-noise ratio with restricted versus amplitude-corrected waveforms 
has already been treated in great detail in 13| (at least for the case of zero spins) so that here 
we can be brief. Let us keep the mass ratio fixed at qm = mi/m2 = 0.1 and let M = mi + 1112 
vary, as in the left panel of Fig. O For low masses the SNR from the restricted waveform 
dominates, but as mass is increased and 2fLso approaches fs from above, the SNR from 
the full waveform catches up. For a distance of 100 Mpc, with the restricted waveforms 
Advanced LIGO can only detect systems with total mass up to ~ 220 Mq and up to double 
that value in EGO; heavier systems have 2fiso < fs and their signals do not enter the 
bandwidth. The inclusion of amplitude corrections effectively doubles the mass reach of 
Advanced LIGO and triples that of EGO. This is because the higher harmonics, having 
higher cut-off frequencies, can still penetrate the detectors' bandwidths for masses where 
the second harmonic does not, and they can lead to sizeable SNRs. 

The right panel of Fig. [2] shows the dependence of SNR on mass ratio for a fixed total 
mass. Here the behavior is qualitatively the same for restricted and amplitude-corrected 
waveforms: In going from = 0.1 to = 1 the SNRs increase by a factor of a few. 



B. Extrinsic parameters 



Let us briefly discuss the extrinsic parameters (tc, ln(l/r), V'c, cos(6), cos(^), V')- 
The last three of these, (cos(6'), 0, ?/;), enter the waveform only through the beam pattern 
functions fl2.10p . and they do so in very similar ways. As a result, these three variables have 
large correlations (in absolute value) with each other and very small correlations with all 
of the other variables. Indeed, using the Fisher matrix based on all of the parameters in 
fl2.22p . typically one has 

^cos{e),4> _ _^cos(e), ip _ _fJ',i> _ 2 (3 1) 

with an accuracy of eight significant digits, and the correlations of these parameters with 
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Figure 3: Errors on coalescence time for varying M and fixed qm = 0.1 (top) and for varying 
and fixed M = 100 (bottom), in Advanced LIGO and EGO. Here and below the distance is 
100 Mpc and angles are as in the previous figure. 



the other eight tend to stay below 10"'^ in absolute value. These angles are unmeasurable 
in practice, with Acos(6') ^ 1, A0 ^ 27r, and Atp ^ 27r. The distance r is also effectively 
unmeasurable, Ar/r ^ 1, as is the phase at coalescence ipc- As mentioned before, Acos(6) 
tends to be reasonably small (~ 0.2 in Advanced LIGO and ~ 0.02 in EGO), but the 
information it carries is not of much interest by itself. 

This leaves the coalescence time tc- As can be seen in Fig. 121 for small (< 50 Mq) masses 
the restricted waveforms give errors of a few tens of milliseconds in Advanced LIGO and 
a few milliseconds in EGO; as total mass is increased, tc quickly becomes unmeasurable in 
practice. By contrast, using the full (2.5, 2.5)PN waveforms, the errors stay below 10 ms for 
M < 250 Mq in Advanced LIGO and below 3 ms for M < 500 M© in EGO.^ For a fixed 
total mass of M = 100 Mq and increasing qm, there is little change in Ate- 



^ At first glance our numerical results for the restricted waveform in botli tliis and tlic following subsections 
may seem at odds with 



26, 



3l|. However, (i) in [26j no spin-related parameters were taken into account, 



which, as shown in 31[, leads to smaller errors; and (ii) in [3l| a somewhat idealized noise power spectral 
density was used which is quite different from the more accurate analytic PSDs now available. Leaving 
out /3 and a we were able to reproduce the results of [26| for both Initial and Advanced LIGO, and those 
of 



3l| were retrieved when using their PSD. 



14 



The overall situation remains qualitatively the same irrespective of mass ratio. Thus, ex- 
cept for the inclination angle, the parameters that give the binary's position and orientation 
with respect to the detector are unmeasurable in practice, both with Advanced LIGO and 
EGO. (Things would of course be different with a network of detectors 32|, but that is out- 
side the scope of this paper.) With amplitude-corrected waveforms the time of coalescence 
is reasonably well determined over a large range of masses, while with the restricted ones tc 
can be measured for stellar mass binaries only. 
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Figure 4: Errors on chirp mass, relative mass difference, and component masses for = 0.1 in 
Advanced LIGO and EGO. 
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Figure 5: Errors on chirp mass, relative mass difference, and component masses for M = 100 
in Advanced LIGO and EGO. 



C. Intrinsic parameters 

We now turn to the errors on the more interesting intrinsic parameters (In(A^), 6, f3, a). 

For the case of fixed g™, = 0.1 and varying M, the errors on Ai are shown in Fig. H] using 
the restricted and (2.5, 2.5)PN amphtude-corrected waveforms in Advanced LIGO and EGO. 
As with tc, using restricted waveforms the errors tend to be small for low total mass but rise 
very steeply with increasing M, becoming effectively unmeasurable as one moves beyond the 
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Figure 6: Errors on the spin-orbit parameter (3 and the spin-spin parameter a with increasing 
total mass, keeping the mass ratio fixed at qm = 0.1. The largest values P and a can take for 
sub-extremal black holes are shown by the horizontal lines. 



mass range of [2Mq, 100 Mq] corresponding to stellar mass binaries. Using the amplitude- 
corrected waveforms they stay below 5% up to M = 250 Mq in Advanced LIGO and below 
1% up to M = 500 in EGO. The corresponding numbers for A6 happen to be quite 
similar when expressed in percent.^ Fig. [5] shows how errors change with varying qm for 
M = 100 Mq fixed. 

What do the errors on Ai and 5 imply for the measurability of individual component 
masses mi, m2? For mi < m2 one has 



The error propagation formula then yields 



Am,- 



rrii 



d In rrii 
d\nM 



V M J 



5-5 



din rrii 
86 



(ASf 



SlnmA f dlnrrii 



dlnM 



d6 



+ 



5(1-52) 



{A5Y + 2 



5(1-52; 



M 



(3.2) 



Y}n{M),, 



(3.3) 



^ Although S is dimensionless, we opted to express its errors in percent. They tend to be small numbers, 
which otherwise would have led to inconvenience when presenting them in tables, as we will do later on. 
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ratio, keeping total mass fixed at M = 100 Mq. Again the largest possible values of (3 and a for 
sub-extremal black holes are shown by the horizontal lines. 



for z = 1,2; in the last term of the RHS the upper sign refers to z = 1. With restricted 
waveforms these tend to be very large 33|]. The situation is dramatically different if am- 
plitude corrections are taken into account, as can be seen in Fig. |H In Advanced LIGO, 
for Qm = 0.1 systems at 100 Mpc, with amplitude-corrected waveforms the relative error 
Ami/mi (resp. Am2/m2) remains smaller than 10% (resp. 15%) for 22 Mq ^ M < 200 M©, 
while with the restricted waveform these errors reach 100% already at M ~ 50 Mq. In EGO, 
with amplitude corrections Arrii/mi (resp. Am2/m2) is smaller than 2% (resp. 4%) for total 
mass up to 500 Mq! 

For a fixed total mass of M = 100 Mq, the relative errors Ami/rrii, i = 1,2, vary signifi- 
cantly with mass ratio Qm (Fig- ED, getting smaller as the system becomes more symmetric. 
For the given total mass, increasing from 0.1 to 1 makes a difference in the errors of 
more than a factor of three in both detectors, bringing the errors down to just a few per- 
cent in Advanced LIGO and a few tenths of a percent in EGO. Note that the errors do 
not perfectly coincide at qm = 1- In the last term in the right hand side of Eq. fl3.3p . the 
correlation a^^^-'^^'^ does not become zero at 5 = 0, and [d\nmi/dS]s=o = —[d lnm2 / dS]s=o 
while [9 In mi /c) In = [dhim2/d\nAi]s=o, which is how the difference in errors comes 

about. This is just an artefact of the quadratic approximation used here to compute errors; 
in the right hand side of (I3.3p all contributions beyond second order are ignored. 
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We end with the spin-related parameters. Again we start with the case of fixed qm but 
varying M. As before the spins are zero so that jS = a = 0. With Advanced LIGO the spin- 
orbit parameter (3 (Fig. [6]) is reasonably well determined with restricted waveforms for lighter 
systems, but it quickly becomes effectively unmeasurable as total mass is increased. The 
amplitude-corrected waveforms do quite well up to M ~ 150 M©. In EGO, the restricted 
waveforms once again lead to reasonably small errors at the low mass end but a quick 
deterioration is seen with increasing M. The amplitude-corrected waveforms lead to Aj3 < 

0. 6 throughout the given mass range. This is to be compared with the largest value |/3| can 
take (~ 9.4) if the binary's components are sub-extremal spinning black holes. Finally, in 
Advanced LIGO, both with restricted and amplitude-corrected waveforms one has Aa ^ 2.5, 

1. e. much larger than the largest value a may be expected to take for sub-extremal black 
holes, which is why we did not include a figure. Also in EGO the spin-spin parameter is 
poorly determined even with amplitude-corrected waveforms. This may well be related to 
the observation by Arun et al. [7| that if the phasing coefficients ^/Ji are formally treated as 
independent parameters, the one at 2PN turns out to be by far the most poorly determined; 
recall that this is the order at which a appears. 

For fixed total mass, the dependence of Aj3 and Act on mass ratio is fairly weak, as 
illustrated by Fig. [71 

Table I: Change in signal-to- noise ratio and improvement of parameter estimation for Advanced 
LIGO and EGO with increasing p in (p, 2.5)PN waveforms. 



AdvLIGO, (5,50)M( 



p 


SNR 


AM/M (%) 


A5 (%) 


Ate (ms) 


AI3 


Aa 





76.5 


13.75 


514.4 


234.3 


51.93 


790.1 


0.5 


85.0 


0.905 


4.498 


6.468 


2.143 


8.529 


1 


74.1 


0.674 


3.662 


5.446 


1.691 


6.908 


1.5 


69.0 


0.463 


4.376 


5.273 


1.124 


6.067 


2 


65.5 


0.458 


3.444 


4.205 


1.144 


5.703 


2.5 


64.0 


0.471 


2.318 


3.822 


1.611 


4.457 


EGO, (10, 100) Mq 


P 


SNR 


AM/M (%) 


A6 (%) 


Ate (ms) 


A/3 


Aa 





461.9 


2.140 


81.75 


74.91 


8.160 


125.4 


0.5 


513.0 


0.145 


0.758 


2.085 


0.347 


1.398 


1 


446.4 


0.107 


0.625 


1.769 


0.272 


1.145 


1.5 


417.7 


0.075 


0.748 


1.753 


0.184 


1.029 


2 


396.8 


0.076 


0.589 


1.405 


0.194 


0.970 


2.5 


387.7 


0.076 


0.401 


1.262 


0.192 


0.755 



D. Dependence on PN order of the amplitudes 

Table [T] shows how errors on (ln(7W), 6, te, (3, a) vary with p for systems with = 0.1 and 
masses for which the SNRs of both the restricted and amplitude-corrected waveforms are 
close to maximum (specifically, (5,5O)M0 in Advanced LIGO and (10, 100) M0 in EGO). 
The restricted, p = 0, waveform contains only the Newtonian harmonic at twice the orbital 
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phase, and the independent parameters characterizing it are (\n{A),\n{Ai),6,tc,ipc, PyCr). 
As explained near the end of subsection IIIBl for p = 0.5 one has two extra harmonics 
and two more parameters become measurable in principle: cos(i) and one of the variables 
{cos{6),(f),ip), although in practice only cos(t) is measurable with any degree of accuracy. 
Each subsequent half-integer increase in p leads to an additional harmonic. However, even 
for p = 2.5 the parameters {cos{6),(j),ijj) are poorly determined and extremely strongly 
correlated, as seen in Eqs. (13.11) . and this is true a fortiori for p < 2.5. 

For those parameters that are measurable, the biggest improvement occurs when going 
from p = to p = 0.5. At p = 2.5 the errors tend to be about a factor of two smaller 
than at p = 0.5. We note that from p = 0.5 onward, all errors generally tend to decrease 
(with minor fluctuations), despite the fact that SNRs decrease monotonically. If we had 
presented errors for fixed SNR, one would have seen an even more pronounced improvement 
in parameter estimation. The great improvement seen in going from p = to p = 0.5 can 
not be accounted for by this being the only case where SNR increases; if we had compared 
these cases for equal SNR then the improvement would still have been very large. Thus, it 
is the extra structure in the additional harmonics that is responsible for the improvement, 
as anticipated in subsection III CI 

For the given system in Advanced LIGO, with the restricted waveform all of the intrinsic 
parameters are unmeasurable (except for Ai, but still with a large error), while with the 
exception of a they are reasonably well determined using the (2.5, 2.5)PN waveform. In the 
case of EGO, note how /? appears to be unmeasurable with the restricted waveform while 
the one at (2.5, 2.5)PN does very well; to a certain extent the same is true of a. 

Table II: Errors for different magnitudes and orientations of the spins with respect to L = (0, 0, 1). 



AdvLIGO, (5,50)M( 



Xl 


X2 


(3 


a 


AM/M (%) 


A5 (%) 


Ate (ms) 


A/3 


Aa 


(0,0,0.9) 


(0,0,0.9) 


8.004 


0.661 


0.449 


2.262 


3.489 


1.314 


4.296 


(0,0,0.9) 


(0,0,-0.9) 


-6.934 


-0.661 


0.494 


2.314 


4.031 


1.214 


4.677 


(0.67,0,0.67) (0,0.67,0.67) 


5.958 


0.366 


0.456 


2.317 


3.581 


1.269 


4.383 


(0,0,1.5) 


(0,0,1.5) 


13.340 


1.836 


0.435 


2.201 


3.350 


1.477 


4.190 


(0,0,1.5) 


(0,0,-1.5) 


-11.557 


-1.836 


0.529 


2.309 


4.253 


1.338 


4.853 


(0.67,0,0.67) 


(0,1.4,1.4) 


12.017 


0.766 


0.438 


2.192 


3.366 


1.428 


4.139 


EGO, 


(10, 100) Mq 
















Xl 


X2 


/9 


a 


AM/M (%) 


A6 (%) 


Ate (ms) 


A/3 


Aa 


(0,0,0.9) 


(0,0,0.9) 


8.004 


0.661 


0.073 


0.387 


1.170 


0.220 


0.732 


(0,0,0.9) 


(0,0,-0.9) 


-6.934 


-0.661 


0.080 


0.391 


1.356 


0.198 


0.786 


(0.67,0,0.67) (0,0.67,0.67) 


5.958 


0.366 


0.074 


0.396 


1.199 


0.211 


0.746 


(0,0,1.1) 


(0,0,1.1) 


9.783 


0.987 


0.072 


0.385 


1.159 


0.227 


0.718 


(0,0,1.1) 


(0,0,-1.1) 


-8.475 


-0.987 


0.080 


0.419 


1.355 


0.214 


0.753 


(0.67,0,0.67) 


(0,1.1,1.1) 


9.527 


0.601 


0.072 


0.385 


1.158 


0.225 


0.714 
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E. Dependence on the magnitudes and orientations of the spins 



So far we have restricted ourselves to systems with zero spins; although (3 and a were 
included as variables in the error analysis and the errors A/9 and Aa were investigated, the 
values of (3 and a were set to zero (i.e., the Fisher matrix was evaluated on the hypersurface 
/3 = (T = 0) until now. In Table HTl we look at the errors for different values of the spin-related 
quantities. The focus is on binary systems where the spins of the component objects are 
close to the Kerr bound, both with < 1 and > 1, z = 1,2. Objects in the latter 
category could be Kerr black holes that exceed extremality, boson stars [s^], or other exotica. 
The spins are chosen at various angles with respect to the orbital angular momentum vector. 
By looking at the results of Table HTl and also comparing with those of Table [T] for p = 2.5, 
one sees that the errors have but a weak dependence on the values of (3 and a. Note that 
AAi/A4 and At^ tend to be somewhat smaller if spins are aligned and larger when they are 
anti-aligned. 

The parameters /5 and a are essentially unmeasurable with restricted waveforms. By 
contrast, with (2.5, 2.5)PN waveforms, in Advanced LIGO and EGO one can measure /3; 
the parameter a remains poorly determined even in EGO. The accuracy with which /3 can 
be measured would already allow for the detection of even relatively small violations of the 
Kerr bound. Apart from spins, the expression (12.131) for (3 depends only on the relative mass 
difference 6, which is rather well determined. Now, if \f3\ is larger than the largest value 
allowed by 6 assuming sub-extremal black holes, then at least one component of the binary 
must have a spin exceeding the Kerr bound. For any given value of S, this "maximum value" 
Pmaxi^) corresponds to the case where Xi ' L = I, i = 1,2, and one has 



113_!^(1_52) 

4 



(3.4) 



We would like to know how well one can measure the difference B{P,6) = \P\ — Pmaxi^)- 
The error propagation formula yields 

{ABf = (A/3)2+ (^^^y(A5)2^2(^^^^ S^'-' 

= {Apr + ^6\A6r T ^5 cr^''A(3A6, (3.5) 
3o 3 

where in the last term, the upper sign is for the case (3 > and the lower one for /3 < 0. 
At P = the two expressions for AB disagree. This discrepancy is analogous to the one 
encountered in subsection IIII CI for Ami/rrii, i = 1,2, ai 6 = 0; again the error is only given 
in the leading order, quadratic approximation. In any case, the last term in (13. 5p tends to 
be small compared to the other terms, including at /? = 0. 

Let us discuss a few examples to probe the extent to which violations of the Kerr bound 
can be measured in Advanced LIGO and EGO. In the above expression for AB, 6 refers to 
the measured value of that parameter; for convenience we set it equal to its physical value. 
Note that in order to claim detection of a Kerr bound violation, it is not necessary for B 
itself to be measurable with great accuracy; all that is needed is the ability to claim B > 
with a reasonable degree of confidence. We take as a benchmark that 2AB be no larger than 
B. 

First let the detector be Advanced LIGO and consider a (5, 50) Mq system (so that 
6 = 9/11) with Xi = X2 = 1.51/. Then one has 2AB = 3.032. This is to be compared with 
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B = 4.447 computed from the actual values of /3 and 6. By the above criterion, a violation 
of the Kerr bound of this size would be detectable. 

Still considering Advanced LIGO and the same component masses, now let Xi = 
(0.67,0,0.67) and X2 = (0,1.4,1.4), setting L = (0,0,1). In this case the first object is 
sub-extremal (||xi|| — 0.95) while the other exceeds the Kerr bound with ||x2|| — 1-98. One 
has 2AB = 2.930 while B = 3.123. Hence also in this case the violation of the Kerr bound 
can be detected. As this example illustrates, if the spins are not aligned with each other 
and L (as will generally be the case), then the Kerr bound will need to be exceeded more 
strongly for the effect to be observable. 

In EGO the parameter (3 is much better determined. Consider a (10, 1OO)M0 system, 
this time with Xi = X2 = I-IL. Here one has 2AB = 0.466 while B = 0.889; hence a 
violation of the Kerr bound at this level can be detected with EGO. 

Finally, take a (10, 100) M© system in EGO with xi = (0.67, 0, 0.67) and X2 = (0, 1.1, 1.1) 
for L = (0, 0, 1), so that one object is sub-extremal with | I — 0.95 while the other exceeds 
the Kerr bound with ||x2|| — 1-56. Then 2AB = 0.462 while B = 0.634, allowing for 
detection of the Kerr bound violation. Once again, with non-aligned spins a stronger excess 
is required for it to be observable. 



F. A note on initial detectors 



Given the above results, it is of interest to similarly compare parameter estimation with 
restricted and amplitude-corrected waveforms in initial detectors, such as Initial LIGO. An 
important caveat is that here the SNRs are expected to be low, in which case the covariance 



matrix formalism may significantly underestimate the errors (25|. We give results for two 
asymmetric systems^ with masses (1.4, 10) Mq and (5, 15) Mq, at a distance of 20 Mpc; 
spins are set to zero. Table UTTl shows errors obtained with the restricted and full waveforms 
for these systems. Again the amplitude-corrected waveform has significant advantages over 
the restricted one. Note that in the case of the lighter system, the full waveform allows for 
determination of (3, at least at the given distance; a is unmeasurable. 



Table III: Parameter estimation with Initial LIGO. Distance is 20 Mpc; spins are zero. 



(1.4, 10) Mo 



Waveform 


SNR 


AM/M (%) 


A6 (%) 


Ate (ms) 


A/3 


Aa 


(0,2.5)PN 


13.49 


1.432 


630.0 


47.47 


12.45 


609.0 


(2.5,2.5)PN 


11.99 


0.262 


29.54 


2.542 


1.223 


29.17 


(5,15)M0 


Waveform 


SNR 


AM/M (%) 


A5 (%) 


Ate (ms) 


A/3 


Act 


(0,2.5)PN 


26.78 


24.72 


4309.1 


101.3 


44.13 


1671.3 


(2.5,2.5)PN 


22.63 


2.598 


19.34 


4.843 


7.590 


20.97 



^ As discussed in subsection IIIDl the Fisher matrix for the restricted waveform becomes ill-conditioned 
when 77ii/m2 — 1 because of the use of (5 as a coordinate, while the Fisher matrix for the full waveform 
becomes singular at mi = 7712 if ln(?7) is used in place of S. For this reason we do not consider the equal 
mass case. 
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IV. SUMMARY AND CONCLUSIONS 



We have compared parameter estimation with restricted versus amphtude-corrected post- 
Newtonian gravitational waveforms, in the stationary phase approximation, for compact 
binary inspiral in Advanced LIGO and a possible third-generation ground-based detector 
called EGO. For the amplitude-corrected waveforms both phasing and amplitudes were 
considered up to 2.5PN order, with inclusion in the phasing of spin-orbit effects at 1.5PN 
and 2PN, and spin-spin effects at 2PN. In the case of restricted waveforms, the parameters 
of interest can be extracted only from the phase. With amplitude-corrected waveforms a 
wealth of information is carried by the spectrum as well, which leads to very significant 
improvements in parameter estimation. Because we mostly dealt with interferometers that 
will detect sources with high signal-to-noise ratios, it was justified to use the covariance 
method for estimating errors and this is what was done throughout. The amplitude-corrected 
waveforms considered here depend on 11 parameters, which makes it difficult to numerically 
invert the Fisher matrix. Fortunately the angles (6', 0, ip) tend to have only small correlations 
with the other parameters, so that it is justified to work with a much better conditioned 
8x8 Fisher matrix in which these angles are disregarded. 

Our main results may be summarized as follows. 

• The parameters that can be measured with reasonable precision using restricted wave- 
forms are 

{ln{A),HM),6,t,,P,a), (4.1) 

although the logarithm of the amplitude, In(^), does not carry much useful informa- 
tion. With restricted waveforms the errors on (14.11) - not counting In(^) - are steep 
functions of total mass, and parameters quickly become unmeasurable as one moves 
towards the higher-mass end of stellar mass binaries. 

Inclusion of amplitude corrections alters the situation dramatically. In that case the 
parameters that can be measured with good accuracy are 

{\n{M),6,t„cos{L),p,a). (4.2) 

Here the errors increase much more slowly as total mass goes up. For sources at 
a distance of 100 Mpc, with amplitude-corrected waveforms parameter estimation 
remains possible for intermediate-mass binaries with total mass up to at least 250 Mq 
in Advanced LIGO and 500 in EGO. Up to these masses, the relative error on 
chirp mass varies the most (by more than an order of magnitude in Advanced LIGO 
and two orders of magnitude in EGO); the errors on the other parameters increase 
by a factor of a few. We also investigated how errors improve with the PN order 
in amplitude; the largest improvement occurs in going from OPN to 0.5PN, although 
apart from fluctuations, errors keep decreasing significantly as the order is increased 
further. At 2.5PN in amplitude, the errors have a relatively small dependence on the 
values of P and cr; A^A/^A and At^ tend to be a little smaller if spins are aligned and 
a little larger when they are anti-aligned. 

• Because our analysis was confined to single detectors, among the intrinsic parameters 
only the time of coalescence can be measured with good accuracy. (The inclination 
angle is also reasonably well determined, but it is not very interesting by itself.) At the 
low-mass end, the full waveform reduces the error on tc by an order of magnitude in 
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Advanced LIGO and about a factor of 5 in EGO compared to the restricted waveform, 
but the differences become much larger as one moves to higher masses. 

• Using the restricted waveforms, component masses tend to be extremely poorly deter- 
mined even for low-mass systems. This is no longer the case with amplitude-corrected 
waveforms: In Advanced LIGO the relative errors on mi, m2 already stay below 15% 
in a significant part of the mass range, while in EGO they are below 4% all the way 
up to a total mass of 500 Mq\ These numbers refer to asymmetric systems; for equal 
component masses the relative errors will be smaller still. The best values seen here 
are as low as a few percent in Advanced LIGO and a few tenths of a percent in EGO. 
All this will have important astrophysical consequences: Unlike what the restricted 
PN approximation suggests, it will be possible to use binary inspiral events to per- 
form black hole population studies, not just for stellar mass black holes but also for 
intermediate-mass black holes that would be the results of successive mergers. 

• The spin-spin parameter a occuring at 2PN in the phasing remains poorly determined 
even when the full waveform is used. With restricted waveforms, the spin-orbit param- 
eter (3 at 1.5PN is undetermined in Advanced LIGO; inclusion of amplitude corrections 
makes it measurable. In EGO f3 is measurable with restricted waveforms for systems 
with low masses, but once again the amplitude-corrected waveforms lead to a very 
significant improvement. On the basis of this we have suggested a way to find out if 
one or both of the component objects violate the Kerr bound. With the full wave- 
forms, the measured value of can be compared with the largest possible value this 
parameter can take given the measured relative mass difference and assuming that the 
binary consisted of two sub-extremal black holes. If for a particular inspiral event, |/3| 
is found to be larger than this maximum value, then at least one component object 
must be more exotic. In this way one would be able to detect fairly small violations 
of the Kerr bound, both with Advanced LIGO and EGO. 

• We also briefly considered parameter estimation in Initial LIGO, showing that there 
too the amplitude-corrected waveform can have a very significant advantage over the 
restricted one. 



Of necessity our analysis was somewhat limited, as the way in which spins enter the 



amplitudes has not been the subject of much attention (see, however, (35|). The drastic 
improvements in the estimation of the other parameters due to sub-dominant harmonics 
suggest that knowledge of spin-related contributions to the amplitudes would hold significant 
advantages. Moreover, in the present work /5 and a were considered constant. As shown 
by Vecchio 3^ and Lang and Hughes 37 1 in the context of restricted waveforms for LISA, 
taking into account the precession-induced modulation due to time-dependent spins also 
leads to qualitative improvements in parameter estimation. 

In both the present paper and its precursor 13| we chose to restrict ourselves to ground- 
based detectors, but it should not be difficult to extend our work to the case of LISA. 



Preliminary studies in this direction have already been performed by Sintes and Vecchio [17 
and Hellings and Moore 



18, 19 



in view of the results presented here it would be of great 
interest to perform a systematic and comprehensive investigation. 

We end with a comment regarding initial detectors. Although in the immediate future it 
would not be computationally feasible to use the full waveforms for inspiral signal searches, 
there would be clear advantages in using them for follow-up studies once a detection has 



24 



been made. As indicated by both our own results and those of Rover et ah 161], this issue 
deserves much more attention than it has received hitherto. 
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Appendix A. Second and third generation interferometric observatories 

The noise power spectral density (PSD) for Advanced LIGO is approximately j26| 



Shif) = So 



-4.14 p. 2 ,111-'- + 



a;-^-^^-5a;~" + lll- 



1 + xV2 



(4.3) 



where x = f / fo, /o = 215 Hz, 5*0 = 10~^^Hz~\ and fs = 20 Hz. The low-frequency cut-off 
is taken to be fg = 20 Hz; below this frequency Sh{f) can be considered infinite for all 
practical purposes. 



Advanced LIGOl 

--■EGO 




f(Hz) 



Figure 8: The strain sensitivities for Advanced LIGO and EGO. 



The basic technical specifications of EGO were outlined in 13|, where an analytic fit of 



its PSD was also introduced. For completeness we repeat the expression for the PSD. Above 
frequencies of a few Hertz it is well-approximated by 



Shif) = So 



X 



1 + bix + b2x'^ + b^x^ + b^x^ + b^x^ + b^x^ 

1 + CiX + C2X^ -|- C^X^ + CiX^ 



(4.4) 
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where again x = f / fo, this time with /o = 200 Hz, and 5*0 = 1.61 x 10 Hz ^. The 
low-frequency cut-off is fs = 10 Hz. One has 

= -4.05, p2 = -0.69, 
ai = 185.62, 02 = 232.56, 
bi = 31.18, 62 = -64.72, 63 = 52.24, 64 = -42.16, 65 = 10.17, = 11.53, 
ci = 13.58, C2 = -36.46, C3 = 18.56, C4 = 27.43. 

(4.5) 

The strain sensitivities for Advanced LIGO and EGO are plotted in Fig. [8l 
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